#!/usr/bin/env python

from math import pi, e, sin, cos
from bernstein import interpolant2D, laplacian

def rhs(x,y):
# Bratu problem-nonlinear
#	y = interpolant2D(n,m,a1,b1,a2,b2,beta,x,y)
#	return e**y
# Poisson equation
        return 2*pi*pi*sin(pi*x)*sin(pi*y)


def lhs(i,j,n,m,x1,x2,y1,y2,x,y):
	return -laplacian(i,j,n,m,x1,x2,y1,y2,x,y)

